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Abstract 

We analyze Regge quantum gravity coupled to SU(2) gauge theory on 4 3 x 

2, 6 3 x 4 and 8 3 x 4 simplicial lattices. It turns out that the window of 

the well-defined phase of the gravity sector where geometrical expectation 

values are stable extends to negative gravitational couplings as well as to 

gauge couplings across the deconfinement phase transition. We study the 

string tension from Polyakov loops, compare with the /3-function of pure gauge 

theory and conclude that a physical limit through scaling is possible. 
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I. INTRODUCTION 



Regge Calculus provides a nonperturbative way for investigations of Euclidean quan- 
tum gravity on a simplicial lattice and offers the possibility to construct a unified theory 
by coupling gauge fields to the skeleton. One remarkable finding was the discovery of an 
"entropy dominated", well-defined phase, where the expectation values with respect to the 
pure Regge-Einstein action are stable The next question addressed was the physical 

relevance of this regime which has been tested by coupling a non-Abelian gauge field to 
gravity PJ. If one assumes that the world without gravity is described by a grand unified 
asymptotically free theory, these numerical studies investigate the relation of the hadronic 
scale to the Planck scale. In particular, it has already been shown that confinement exists 
in the coupled system ||. In this work, we perform a nonperturbative analysis of the phase 
diagram of Regge quantum gravity coupled to SU(2) gauge fields on several lattice sizes in 
four spacetime dimensions. The stability and boundary of the well-defined phase is investi- 
gated on lattices of sizes up to 6 3 x 4, considerably extending the previously available data. 
Within this phase we find that the confinement-deconfinement transition of conventional 
lattice gauge theory is still present. We extract values for the string tension and gain some 
evidence from the /^-function that the window of the well-defined phase extends to large (3 
values corresponding to small lattice spacings a in the order of the Planck scale. 



II. ENTROPY DOMINATED PHASE 

In Regge Calculus the edge lengths are considered to be the dynamical degrees of free- 
dome of the discretized spacetime manifold. In d = 4 dimensions the geometry is Euclidean 
inside of a rf-simplex and the curvature is concentrated at the d — 2-subset of the lattice, the 
triangles. Quantization proceeds via the path integral although the choice of the gravita- 
tional measure is an unresolved issue. Investigations on regular triangulations do not favor 
any of them ||], but for irregular triangulations a preference for scale- invariant measures was 
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found 0. It may be that measures are divided into universality classes such that identical 
physics is obtained within one class. 

In our simulations we chose a hypercubic triangulation with N% x N t vertices and the 
scale-invariant measure: 

m 2 }] = I[^r, (2-1) 

I 1 

where / is used to denote the link label as well as its length. The system of SU(2) gauge 
fields coupled to quantum gravity has the action || 

S = 2m 2 p £ A t a t - | £ W t Re[Tr(l - U t )\ , (2.2) 
t * t 

with the first term being the Regge-Einstein part composed of the bare Planck mass m p and 
A t , a t the area and the deficit angle of the triangle t. The addition of the SU(2) gauge term is 
straightforward and follows ordinary lattice gauge theory by assigning SU(2) matrices to the 
links. The elementary plaquettes become triangles on the Regge skeleton. f3 corresponds to 
the inverse gauge coupling and the weights W t = const x t^at, with a 4- volume V t assigned 
to each triangle, describe the coupling of gravity to the gauge field. U t is the ordered product 
of SU(2) matrices around t. In contrast to the flat lattice, the unit matrix in the action is 
important because the weight factors are dynamical. They are constructed such that the 
correct continuum limit is ensured in the limit of vanishing lattice spacing ||. 

We present Monte Carlo (MC) results concerning the boundary and stability of the well- 
defined phase. For this purpose large statistics simulations were performed on 4 3 x 2 and 
6 3 x 4 lattices, covering a variety of (m^, /3) values. The accumulated statistics is summarized 
in Tables | and fl]. 

The stability of the coupled system was analyzed from MC-time histories of the Regge 
action. Examples are presented in Fig. [[]. Long runs on the larger lattice show that it is 
difficult to decide whether the well-defined phase is stable or just metastable. In the latter 
case an additional curvature term of higher order |J may stabilize the system. 

The deconfinement transition was studied from the behavior of the Polyakov loops 
P = ±Tr(UxU 2 ...U Nt ). For 4 3 x 2 and 6 3 x 4 lattices Fig. @ shows MC-time histories of 
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P and / in the confined and deconfined phase for gravitational couplings in the well-defined 
phase. The link lengths are largely independent of fluctuations of the order parameter. 
Notable is the long equilibration time in the deconfinement phase. 

Extracted phase diagrams of the gauge-gravity system are displayed in Fig. |3|. The dotted 
lines are to guide the eyes and rely on the depicted stable versus unstable data points. From 
the 4 3 x 2 lattices we have numerical evidence that the stable phase extends to (3 = 3.0 
(Table H). The 6 3 x 4 lattices indicate a glitch in the well-defined-ill-defined boundary when 
passing from confinement to deconfinement. However, the present runs do not decide the 
question conclusively. It may be accidental that at (3 — 1.6 the = 0.025 and = 0.023 
data did not run away. For f3 = 1.55 tunneling into the ill-defined phase happened for 
these nip values only after more than 100k sweeps. Our systems exhibit a small shift of the 
deconfinement phase transition from (3 c (N t = 2) = 1.525 to f3 c (N t = 4) = 1.575 with error 
bars less than 0.025. 



III. STRING TENSION 

The Polyakov loop P{R) in the short extent L t of the lattice describes the propagation 
of a static quark. We introduce a quark source and a sink separated by a distance R and 
calculate the correlations of the Polyakov loops at these points. As in conventional SU(2) 
lattice gauge theory at finite temperature T [0)0] > we extract from the correlation function 

(P(0)P\R)) = eM~V(R)}, (3.1) 

the quantity V(R) corresponding to the potential between the static quark-antiquark pair. 
In the confinement phase V(R) should grow linearly for large R due to an infinite free energy 
of isolated quarks: 

V c (R) = ^- + aR + C, (3.2) 
where a is the Coulomb parameter, a the string tension, and C a constant. 
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The correct distance between two points should be measured using geodesic distances. 
We take the distance R between the source and sink to be equal to the index distance 
along the main axes of the skeleton. This seems a reasonable approximation in the well- 
defined phase with small curvature fluctuations. Using scalar field propagation one may 
calculate corrections which are expected to be small for our purposes. 

To extract a reliable value for the string tension, we simulated the system on an 8 3 x 4 
lattice. Our data rely on 30000 measurements after equilibration for the coupled system as 
well as for the pure gauge system on a flat simplicial lattice. 

Figure f|(a) presents the data points for the confinement potentials for several gauge 
couplings in the presence of gravity with ml = 0.005 in the "entropy dominated" phase, 
while Fig. £|(b) depicts the situation with gravity switched off. The dotted lines correspond 
to fits to the correlations in Eq. ( |3.1|) according to the potential of Eq. ( |3.2|) with the Coulomb 
parameter fixed to a = ^ and a mirror term included. 

Since we have extracted string tensions for several f3 values, we are in a position to study 
its scaling behavior. For pure gauge theory /3-functions are derived in the literature also for 
a simplicial lattice |T3J. We fit our string-tension data to the function 



102/121 

Afl a t V 11 ) ^ V 11 



a phys 102/121 /_ 6ttV5^ 

. n „ „ exp _ _ 



a flat = ^ exp -—^ (3.3) 



and obtain a value for Afl at = 0.0102(l) v /a p hys- This is in good agreement with an analysis 
of Wilson-loop ratios in the T = case, yielding A flat = 0.008(l) v /cr p hy S . For the fluctuating 
lattice to our knowledge a /^-function is not worked out, we are aware only of a recent study 



for the pure gravity case within dynamical triangulation [fL4| . Thus, we used as a starting 
point the above SU(2) function Eq. ( |3.3|) . In Fig. [| we compare the string-tension scaling 
for the flat and the fluctuating lattices. 

To have both systems on the same scale, we rescaled the inverse gauge coupling 
[3 — > (W t )f3 of the fluctuating system. We find A grav = 0.0090(l) v /cr p hy S which is very similar 
to the flat case. Assuming that the pure-gauge /3-function is a reasonable approximation for 
the full renormalization relation, our data show that the scaling window opens already at the 



(3 values considered, similar to the pure SU(2) case [|13|]. As a consequence, the continuum 
limit could be performed along (3 —* oo, eventually stopping at a finite value corresponding 
to a lattice spacing equal to the Planck length. The existence of a corridor to large /3's is 
indicated in the (m p , (3) diagrams of Fig. |^. 



IV. CONCLUSION 

The boundary between the well-defined and the ill-defined phase of Regge quantum grav- 
ity coupled to SU(2) gauge theory was studied with the largest so far available statistics on 
6 3 x 4 lattices. These lattices are already very CPU time intensive, because the gravitational 
dynamics is very slow. Altogether, evidence was gained that the well-defined phase is stable 
with the increase from 4 3 x 2 to 6 3 x 4. 

Within the well-defined phase we find that the confinement mechanism from the non- 
Abelian gauge fields is not spoiled by quantum gravitational effects. This is not trivial, 
because it is was not clear how quantum gravity affects a gauge theory. In the confined 
phase we observed a potential linearly rising with R. Extracting string tension values for 
both the coupled system and for the pure gauge theory on a simplicial lattice without gravity 
we found a very similar scaling behavior. It indicates that gravity effects do not destroy the 
physics of conventional asymptotically free field theories, even if the gravity-gauge coupling 
is large as in our situation. This gives hope that the physical limit can be approached 
through scaling from the investigated region towards the Planck length. Additionally, to 
reproduce Newton's law one should demonstrate a diverging gravitational correlation length 
of graviton propagators when approaching the boundary of the well-defined phase. 
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FIG. 1. MC-time histories of the Regge action (A t St) for 4 3 x 2 (upper) and 6 3 x 4 (lower) 
lattices, with one time unit representing lk sweeps. The Planck masses used in the simulations 
are indicated in brackets beside the curves. The plots correspond to simulations at the (3 values 
given in the titles. 
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FIG. 2. Time histories of Polyakov loops P and link lengths / on 4 3 x 2 (upper) and 6 3 x 4 
(lower) lattices, with one time unit representing lk sweeps. For reasons of comparison link lengths 
are normalized with respect to the maximum of all four runs. In the lower figures P is multiplied 
by a factor of two. The (m^, (3) values are given in the titles. 
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FIG. 3. Phase diagrams extracted from 4 3 x 2 (a) and 6 3 x 4 (b) lattices. The bold vertical 
lines indicate the confinement-deconnnement transition. 
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FIG. 4. Static quark potentials on an 8 3 x 4 fluctuating lattice with nip = 0.005 (a) and for a 
flat simplicial lattice with gravity switched off (b). The curves are fits to a confined potential with 
periodicity taken into account. Error bars arise from a jackknife procedure. 
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FIG. 5. Scaling of the string tension a for the flat and the fluctuating lattice (nip = 0.005). 
a was extracted from fits to Eq. (3.2) using a = it/ 12. 
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TABLES 



TABLE I. Statistics for the 4 3 x 2 lattices in units of lk sweeps. 
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TABLE II. Statistics for the 6 3 x 4 lattices in units of lk sweeps. 
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